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By R. F. Hoelker and B. P. Winston 
Electronics Research Center 


SUMMARY 

In two concurrent series of graphs, a class of orbits in 
the Earth-Moon (E-M) field and a class of Kepler orbits in rota- 
ting coordinates are depicted and compared. 

A general discussion of characteristics of rotating Kepler 
orbits is included. 

The model used for the E-M orbits is that of the restricted 
problem of three bodies. The orbits of the class depicted origi- 
nate at the E-M line, half of the E-M distance beyond the moon 
with velocity orthogonal to the E-M line within the E-M plane. 


INTRODUCTION 

This report is the first of a series intended to depict 
classes of orbits of the restricted problem of three bodies and 
to compare these orbits with orbits of the Kepler problem, the 
latter of which refer to a rotating coordinate system. 

The purpose of these publications is twofold: Since the 

model of the restricted problem with suitable mass-ratio may 
serve as a good approximation to the Earth-Moon (E-M) model, the 
results will be of assistance to the mission designer concerned 
with the selection of E-M trajectories. In accordance with this, 
the terms "Earth-Moon orbits" and "orbits of the restricted three 
body problem" will be used interchangeably in this report. 

The second purpose is tutorial — that of helping the 
scholar in acquiring familiarity with the characteristics of 
orbits of the restricted three-body problem, and, by the con- 
current showing of comparable Kepler orbits, contributing to the 
understanding of many peculiarities encountered in the restricted 
problem. 

The class of E-M orbits discussed in this report is deter- 
mined by the common initial position on the axis of the masses, 
one half E-M distance beyond the moon and the initial velocity 
being orthogonal to the mass-axis, in the plane of motion of the 
masses. The range in which velocity variations are made is 
essentially that range for which the Kepler orbits are elliptical 



Since Kepler orbits , as viewed from a rotating reference 
system, have obtained rather limited promulgation, an introduc- 
tory section on these orbits is included as part of this report. 

The comparative study of the orbit classes is arranged in 
two sections, the first offering a synoptical viewing of orbit 
arrays with each orbit shown for a limited length of time, and 
the second section tracing and comparing single orbits over ex- 
tended histories. 

The appendix compiles information related to the mathe- 
matical problem description and the methods of regularization 
and integration as applied in the computational program that was 
used in the preparation of the material presented. 


KEPLER ORBITS IN ROTATING COORDINATES 

The model of the Kepler problem is arrived at by making the 
following transition from the normalized model of the restricted 
problem of three bodies. The smaller of the two attracting masses, 
i.e., the mass p, is reduced to zero while the sum of the two 
masses is maintained at unity. In this process both the distance 
between the two attracting masses and their rate of circular revo- 
lution about their common mass-center are maintained at unity. 

Two Cartesian coordinate systems, the inertial (X, Y) and the 
rotating coordinate system (X R , Y R ) , are used in this report. 

The latter rotates with the masses, the X R -axis placed through 
the two attracting centers and the Y R -axis placed through the 
common center in the plane of their revolution. The two coordi- 
nate systems are coinciding at initial time. 

After the transition from the model of the restricted three- 
body problem to that of the Kepler problem has been made, the 
origin of both systems is at the only remaining mass. 

The limit point of the location of the second (vanishing) 
mass is at the X R -axis at distance unity. This point may be 
coined the "empty mass point" and marked "L" for later reference, 
where its significance will be discussed. As to the orientation 
of orbits in the Kepler problem, the convenience of having the 
apses of the orbits on the abscissas of the two systems at initial 
time is utilized. In this way, the objective of comparing these 
orbits to those of the three-body problem is met most easily. For 
the explanations to follow, attention is drawn to the different 
definitions that some terms will have in the two systems of refer- 
ence used here, i.e., the inertial (or space-fixed or sidereal) 
and the rotating (or synodic) coordinate systems. Frequent use 
of the identification of orbits by their mean-angular motion "n," 
as defined with reference to the inertial system, will be made. 
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The reason for this lies in the convenient interpretation this 
term finds in the rotating system. This will now be developed. 

A condition for Kepler orbits to be periodic or returning 
to their initial state coordinates in finite time in the rotating 
reference system is the commensurability of their inertial period 
Pj with the period 2ir of the rotating system in the inertial 
frame. Making use of the relationship n = Pj/2tt results in the 
requirement for n to be rational for orbits to be periodic in 
the synodic system. 

If one assumes n = k/i , with k and l being relatively prime 
integers, the following interpretation of n is obvious. During 
one synodic period of the orbit with characteristic n = k /£,, the 
coordinate system rotates i-times, while an orbiter moving along 
the given orbit revolves k-times along this orbit, with both 
motions understood to be in reference to the inertial frame. 

The term k is signed positive or negative in agreement with 
the sense of motion along the orbit as seen in the inertial frame. 

Since information on synodic Kepler orbits, especially 
their illustration, is rather sparse in the literature, this re- 
port will devote space to a systematic discussion of them. (See 
also Ref. 1.) 

Figure 1 shows the inertial representation of a group of 
Kepler orbits all of which have mean angular motions of n = 3/4. 

A fixed n-value for a group of Kepler orbits means as much as a 


Y -AXIS 



Figure 1.- Kepler orbits with n = 3/4 in inertial coordinates 
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fixed major axis, fixed period, and fixed energy for this group. 
The notion "family" is often applied to the totality of orbits 
of a fixed n. Note that the extremes of a family are the rec- 
tilinear flights along the X-axis between the point X = 2a and 
the point-mass at the origin. 

The two subsequent illustrations (Figures 2 and 3) depict 
these orbits and several others of this family as they are viewed 
from the synodic coordinate system. Of the original geometric 
shapes there is only one retained, the circular orbit marked 
"G" on Figure 1. Inspecting now the series of orbits, starting 
at the circular one and moving in either direction, reveals a 
pattern development that is formed sequentially by indentations, 
cusps, and then loops of increasing magnitude. Each loop develops 
on its own an increasing indentation where it is nearest to the 
mass. The series terminates at both ends with orbits that show 






Figure 2.- Kepler orbits with n = 3/4 
in rotating coordinates; Part 1 



(COLLISION) 
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a cusp-like geometry at the mass (orbits "A" and "N"). These 
latter cases are the collision cases, corresponding to the recti- 
linear cases in the inertial reference. (The term "cusp" refers 
to the geometry only. The infinite velocity at collision actually 
prohibits the use of the term here and so, subsequently, the term 
"cusp" will be applied only where both velocity components — 
with respect to the rotating coordinate system — vanish.) 


Just as the rectilinear cases in the inertial reference can 
be associated to the positive-n family as well as to the negative- 
n family, so can the last discussed collision cases in the synodic 
reference be considered to terminate the positive-n family and to 
commence the family of the negative n-value. Synodic orbits of 
the negative-n family are represented on the subsequent illustra- 
tions (Figures 4 and 5). The first orbit ("0") and the last one 
("X") are repetitions of orbits "N" and "A," respectively. 




X R = 2.3; Y R = -2.51 


ORBIT Q 




X R = 2.0; Yr =-2.42 


X R = 1.75; Yr = -2.31 



X R = 1.50; Yr = - 2. 21 


Figure 4.- Kepler orbits with n = —3/4 
in rotating coordinates; Part I 


Figure 5.- Kepler orbits with n = —3/4 
in rotating coordinates; Part II 



(CIRCULAR) 
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Tracing the pattern development from these two extremes in direc- 
tion to the center of this series, allows one to recognize that 
loops about the mass are formed which increase in size as the 
center of the series is approached. The contraction of the outer 
contours of the orbits which can be observed to occur with the 
expansion of the inner loops leads finally to the coinciding of 
the outer and inner traces to form a circle (orbit "T" ) . (Note 
that, for ease of reading, some of the orbits are traced for a 
third of the period only, the two other thirds being identical to 
this, but rotated by 120 and 240 degrees.) With the attainment 
of this circle, the full cycle, consisting of two families, is 
closed. 


It is worth mentioning that the two circular orbits, "G" 
and "T," evolve from neighboring patterns that are topologically 
quite different from each other. While the circular orbit "G," 
which is of the positive family, matures through the process of 
contraction of loops and smoothing of dents, the orbit "T" of 
the negative family is formed through the convergence of a number 
of loops about the mass. Since here in the negative n-range, 
the mass center is circumscribed seven times by all orbits of 
the family, this circular orbit is to be run through seven times 
if consistency of its orbital period with that of its neighboring 
orbits of the family is to be maintained. The orbit G, in con- 
trast, fits into its neighborhood if traced through only once. 

To further enhance the understanding of, and familiarity 
with, Kepler orbits in the synodic system, five more families 
are depicted on the five subsequent illustration^' (Figures 6-10) . 



Figure 6.- Kepler orbits periodic in the 
rotating system: positive; n= 4 


Figure 7.- Kepler orbits periodic in the 
rotating system: positive; n - 2 
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These concern the n-values 4, 2, 3/2, 4/3, and n = 1, all of 
positive sign only. (Note that, in most cases, several orbits 
have been drawn on one grid. Also, for the family n = 4/3, the 
limiting orbits, i.e., the collision orbits, are traced to one 
loop only while four of these make up the orbit.) 

Characteristic for orbits of positive n-values that are 
larger than unity is a geometry marked by loops extending "out- 
ward" in contrast to the loops (rather lobes) of the (n < 1) 
cases which are directed "inward." 

A case of uniqueness and, thus, of special interest is the 
family of n = 1, depicted on Figure 10. Since for this family 
the mean angular motion is the same as that of the rotating 
system, the orbits do not enclose the origin. Thus there also is 
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no circular orbit in the family. Rather the circular orbit in 
the inertial reference has its synodic image in a stationary 
point at distance unity (X R = 1.0), i.e., the point marked "L" 
before. All orbits of the family envelop this point. Also, 
orbits of this family do not intersect one another. (For a 
discussion of these orbits, see ref. 2.) 

Families of orbits have been discussed with respect to 
their internal development from member to member. An attempt 
will now be made to compare families with each other and thereby 
identify the characteristic features of families. The larger 
generality of statements then calls for omission of particulars 
which, in this case, is accomplished by selecting from each 
family only one representative member. The selection of each of 
these members is based on its capability of demonstrating the 
typical features of its family, as the number of loops, their 
structuring and sequencing, the number of times the orbit is 
going around the mass in a period, and so forth. Of the families, 
15 have been chosen. They are shown as Figures 11-14. Again, 
positive-n families only are shown since the structure of negative- 
n families is rather easily understood. The variation of n among 
the illustrated cases is broad enough that the behavior of any 
other n-case can be found by analogy. These illustrations are 
evaluated in the form of a list of characteristics which follows. 




Figure 12.- Kepler orbit characterization 
(Part II) 


Figure 11.- Kepler orbit characterization 
(Part I) 
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Characteristics 


The following statements hold for all positive and negative 
n-families, except a few cases listed at the end of this section. 
The proof of the characteristics is foregone; the accompanying 
illustrations may serve to make verification and to gain under- 
standing . 

For Kepler orbits of mean angular motion n = k/Z with k 
and l being relatively prime integers, k ^ 0 and Z > 0 , the 
characteristics of orbital motion, as observed in a synodic (i.e., 
rotating) coordinate system, are the following: 


(a) The synodic period is 2nZ. 

(b) The mean-angular motion with respect to the synodic 
reference system is (k - Z)/Z. A negative value means 
retrograde motion. 

(c) During one synodic period an orbiter moving along the 
orbit goes around the mass (k - Z ) -times as observed 
in the synodic reference system. 
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(d) During one synodic period, the orbit develops jkj 
identical patterns (loops, lobes, cusps, etc.) 

(e) Two patterns that succeed each other in time are spaced 
from each other by a central angle of 2ir(k - £)/|k|. 

(f) Families of negative n will not have a cusping orbit. 

Qualifications are in order for the application of these 
rules in the following cases : 

1. Circular Orbits 


The synodic period according to rule (a) need not be the 
smallest one for circular orbits. If |k - £ | > 1, an earlier 
period is calculated by 2 tt£/ | k - £|. However, in this case, 
consistency with the period of the family will not be maintained. 
A corresponding remark holds, of course, for rule (c) . 

2 . Parabolic Orbits 

Though parabolic orbits can be thought of as associated 
with n = 0 and will be marked correspondingly, their character- 
istics are found by approaching them through a sequence of small- 
n-families as, e.g., the sequence l/£ with ever-increasing £. In 
this way, one finds that the period is infinite and the mean 
angular motion in the synodic system is negative unity. This 
holds for the positive as well as the negative parabolic family. 

3 . Libration Orbits and Libration Point L 


The family of Kepler libration orbits (n = +1) possesses 
neither a circular nor cusping orbits. These are represented in 
degenerate form by the Kepler libration point L. 

The assignment of a period to the stationary point L, 
according to rule (a) , would only serve the purpose of preserving 
consistency throughout the family. 


THE APSIDAL VELOCITY DIAGRAM 

Since the approach to the three-body problem in the suc- 
ceeding part of this report will not be in the form of a search 
for periodic orbits that make up a particular family, but will be 
rather in the form of a study of orbits that all start at a fixed 
origin — and in a comparison of these orbits with Kepler orbits 
that also start all at one fixed origin — a diagram which relates 
initial conditions to the n-value of an orbit is needed. Such a 
diagram will now be introduced (Figure 15) . 
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Figure 15.- Apsidal velocities in both systems for families of orbits with fixed 
mean motion **n M versus apsidal position on X-axis or Xp —axis 


A Kepler orbit is determined by its initial conditions of 
position and velocity. Let the Y-value and X-value be zero at 
initial time for all orbits under investigation. Then any orbit 
can be represented as a point in the system, . the abscissa of which 
is X and the ordinate of which is Y. Since Y is the velocity at 
the apsidal point X, the representation is called the apsidal 
velocity diagram. If orbits of identical mean-angular motion n 
are connected by lines in this diagram, the desired relationship 
is furnished. 

Such a representation is given in the left-hand graph of 
Figure 15. Families of n = 2 , 1, and 1/2 are shown. Since here 
the coordinates are inertial, the curves are symmetric with 
respect to the X-axis and the sign of the Y-value is equal to the 
sign of the n-value. 

If the fact is borne in mind that orbits along a constant-n 
curve also have constant major axis 2a, it will be recognized 
that the intersections of the n-curves with the X-axis give the 
respective rectilinear orbits, starting at 2a, colliding with the 
mass and returning to the initial point. Also, it will be recog- 
nized that the heavy dot marks on each curve, at an X-value half- 
way toward the maximum of the curve, reprsent the positive and 
negative circular orbits of the family. 

On the right-hand graph of Figure 15 the same three n- 
families are shown, now plotted with respect to the rotating 
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coordinates Xr and Yr. Since at initial time X = X R , the 
abscissa-values are unchanged. The ordinate-values transform 
according to Y r - Y - ojX, with w being the angular motion of the 
synodic system of value unity. 

Note here that now the major part of each family is initi- 
ated with negative velocity. Also all collision cases (lying on 
the dashed straight line) are "in the negative," and so are both 
circular orbits for initial positions larger than unity. Observe 
also that, on the lower branches of some n-curves , X R is double- 
valued over some regions of Yr. 

While the two graphs of Figure 15 serve primarily to intro- 
duce the concept of the apsidal velocity diagram, the diagram 
shown in Figure 16, by preference, will be used as the working 
diagram in the following discussion. Its more detailed structure 
shows the field of n-curves wedged in between the two zero n- 
curves. These limits represent all orbits of parabolic character. 
The highest n-value shown here is n = 2, though, of course, the 
field continues to the left of this curve. 



Figure 16.- Apsidal velocities in the rotating system for orbit families with fixed mean 
motion versus apsidal position on Xp —axis 

Principal interest will be directed to the intersection of 
the n-field with the vertical line given by Xr = 1.5 since all 
the Kepler orbits depicted and discussed in this report will be 
initiated at this distance from the mass. Also, the orbits of 
the restricted problem of three bodies, discussed in this report, 
will originate at a distance of 1.5 from the heavier body. 
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The discussion of the Kepler series in the next section 
will commence with an orbit near the lower parabolic limit, will 
carry through the elliptical range, and terminate when orbits of 
small positive apsidal velocity are reached. As can be seen in 
this diagram, the n-families to be crossed range from n = 0 
(negative) to the value near n = -3/2 whereupon the collision 
trajectory will be hit. Subsequently the n-values are positive 
and reverse in order from n = +3/2 down to n = 0 (positive) . In 
the series, two circular families (dashed curve) will be crossed 
and two orbits with cusps — one of which will be the orbit with 
zero initial velocity — will be found. 

While this apsidal velocity diagram will serve as a guide 
for the Kepler orbits as the study progresses, so the Kepler 
orbit series as a whole will serve as an approximation to the 
series of three-body orbits started at a corresponding initial 
point. Both series will be discussed simultaneously. 


COMPARATIVE GRAPHICAL DESCRIPTION OF A SERIES 
OF ORBITS OF THE RESTRICTED PROBLEM OF THREE BODIES 
AND A SERIES OF KEPLER ORBITS IN ROTATING COORDINATES 


Preparatory Remarks 

The model of the restricted problem of three bodies, as 
employed here, is that of two point-masses (1 - y and y) , at 
unity distance from each other, revolving with unity angular 
motion oo about their common mass center. The mass parameter y 
is chosen as 1/80. A rotating Cartesian coordinate system Xr 
and Yr is defined such that the origin is at the mass center and 
the masses on the abscissa at Xr = -0.0125 and +0.9875 for the 
larger and smaller masses, respectively. For convenience, the 
masses will be referred to frequently as Earth and Moon, or E-M. 

The objective of this report is to bring into comparison 
a series of orbits of the restricted three-body problem with a 
series of Kepler orbits. The three-body series is defined by 
the initial conditions 


X R = 1.4875; Y R = 0 ; X R = 0 ; Y R varied, 


the differentiation being understood with respect to the normal- 
ized time variable. The Kepler series collated to this series 
is defined by Xr = 1.5000; Yr = 0; Xr = 0; Yr varied. The initial 
points are chosen such that their distances to the heavy mass (or 
only mass, in the Kepler case) are identical in the two series. 
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There is no fixed time length for which orbits are studied; 
rather the time is dictated by the case history itself. On the 
average, however, the time length lies between 2v and 4 it, i.e., 
the time of one and two lunar revolutions about the Earth. 

It will not be a matter of primary concern in this study to 
search out periodic orbits. However, in most cases where, in the 
course of progression, cases of periodic orbits are encountered, 
they are isolated and illustratively documented. In other cases, 
it was considered satisfactory to point out the existence of 
periodic cases by showing sequences of orbits that clearly are 
positioned at both sides of a periodic orbit. Where a periodic 
orbit is pictured, it is usually traced to its first orthogonal 
crossing of the X R -axis. This crossing of the axis is at half 
the orbit's period, the second, reflected half of the orbit being 
omitted in the interest of greater clarity. 

The comparison of the two series will be made twice. The 
first time, an overall view of the orbit series of one problem 
will be set in correspondence to an overall view of the counter- 
series. This view will be concerned with short time histories of 
the orbits only and therefore will allow the composite exposition 
of a fair number of orbits on a single graph. 

In the second comparison, individual orbits of the restric- 
ted problem will be depicted one by one, in a sequence that will 
be dense enough to relate pictorially all events of relevance as 
far as they occur in the time interval of 4n after the initial 
time. These will then be set into relation to Kepler orbits as 
they fit best. 

Subsequently, the three-body orbits will be consistently 
shown at the upper half of the page and the Kepler orbits on the 
lower half of the page. 


Comparison of Short-Term Synopses 

The near-time behavior of the orbits of the restricted 
three body problem, starting at X R = 1.4875 in axis-orthogonal 
direction, is aptly characterized in Figures 17, 19, and 21. 

Their Kepler counterparts are shown as Figures 18, 20, and 22. 

The upper series (Figure 17) is started with the orbit of 
initial velocity Yr = -2.47, a value not far from the escape cases 
which would be represented by orbits spiralling out to infinity. 
Initial velocities on this first graph are increasing to a value 
of Y r = -1.5, associated to the innermost orbit. Intersections 
of the orbits with the Xg-axis move successively closer to the 
Earth. The last orbit of this group is an Earth-collision orbit. 
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PART I : INITIAL VELOCITIES f R 
BETWEEN -2.47 AND -1.50 



PART H : INITIAL VELOCITIES t R 
BETWEEN -1.50 AND -1.06 


Figure 17.- Earth-Moon field: orbit synopsis, Part I and II 

The next group of Figure 17, which repeats the last colli- 
sion orbit and shows orbits of initial velocities up to Y R = -1.06, 
demonstrates how the orbits "squeeze-in" between the two masses 
to intersect the X^-axis. The second intersection of the orbits 
with the XR-axis, however, moves further out toward the point 
where the orbits start. Studying then the corresponding graphs 
of Figure 18, one observes a close similarity of the Kepler orbits 
with those of the three-body problem. 

The obvious difference of behavior develops when orbits of 
the upper series collide with the Moon, which event occurs once 
in each of the two groups. It will be shown in the next section, 
that there is a neighborhood of orbits to these collision orbits 
that will show the influence of the moon on their flight course. 




PART I: INITIAL VELOCITIES f,, 

BETWEEN -2.47 AND -I.S0 


PART n: INITIAL VELOCITIES * R 

BETWEEN -1.50 AND -1.00 


Figure 18.- Central force field: orbit synopsis. Part I and II 
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PART m ■■ INITIAL VELOCITIES 
BETWEEN -0.980 AND -0.849 



PART 12 : INITIAL VELOCITIES 9r 
BETWEEN -0.8488 AND -0.7600 


Figure 19.- Earth-Moon field: orbit synopsis, Part III and IV 


Progressing now to Figure 19, it can be seen that from 
orbit to orbit the loops around the Moon contract to a smaller 
size until lunar impact is encountered with the initial velocity 
of -0.849. The orbits that follow this collision orbit in the 
sequence of the series, as shown in the second graph of Figure 19, 
circumnavigate the Moon in a direction opposite to that flown by 


the orbits on the left-hand graph 
the region ahead of the Moon, 
they return to the original di- 
rection of mean angular motion. 
With increasing initial velo- 
city, the loops of the orbits 
shrink and finally degenerate 
into a cusp. 

If there is interest in 
large-scale characterization of 
the orbits, e.g., a character- 
ization of entire groups among 
each other, the two groups of 
Figure 19 are both concerned 
with "near-Moon" orbits, while 
this characterization is not 
applicable to the two first 
groups discussed. As will be 
seen in the next section, the 
"near-Moon" groups have among 
their orbits a number of peri- 
odic ones that stay permanently 
in the neighborhood of the Moon. 


after a larger "excursion" into 



INITIAL VELOCITIES Yr 
BETWEEN -0.962 AND -0.758 


Figure 20.- Central force field: orbit synopsis, 
Part III 
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PART 2: INITIAL VELOCITIES Yr BETWEEN 
0.738 AND +0.30 


Figure 21.- Earth-Moon field: orbit synopsis, 
(concluded) 



INITIAL VELOCITIES Y R BETWEEN 
-0.738 AND +0.30 


Figure 22.- Central force field: orbit synopsis, 
Part IV (concluded) 


The group of Kepler or- 
bits that forms the counter- 
part to the "near-Moon"-groups 
of the three-body problem is 
depicted in Figure 20. The 
development is rather clear- 
cut, with loops decreasing in 
size and simultaneously moving 
in retrograde direction. At 
the end of this development, 
there is again the orbit with 
cusps, with nearly the same 
initial velocity (-0.758) as 
was used for generating the 
cusping orbit of the three- 
body series. The Kepler series 
development is simpler than 
that of the three-body series 
since the disturbance by the 
Moon is absent. 

Figure 21 then depicts 
the terminal phase of the orbit 
development of the Earth-Moon 
field. The group starts at 
the innermost orbit, which lies 
in the development near to the 
cusping orbit of the former 
group. With the initial velo- 
city increasing, the orbits 
assume patterns that enlarge 
radially from one orbit to 
the next. The general develop- 
ment, one notices, is the re- 
verse of that of the very first 
group of the sequence (i.e., 
that of Figure 17) , which 
means that finally orbits are 
obtained which are spiralling 
outward at a rate that grows 
with the initial velocity. 

The last orbit shown (with 
Tr = +0.3) is well in the 
range of escaping orbits. As 
Figure 22 shows, this develop- 
ment is not different from 
that of the Kepler series . 
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Comparison of Individual Orbits 
over Longer Periods 

Both series will now be considered once more, this time 
in a denser orbit progression and with attention to longer 
orbital histories. 


The density of selection of orbits in the Kepler series 
will be dictated by the objective of presenting a coherent 
development of the family-structure in the same sense as men- 
tioned in an earlier part of this report. 


For the series of the restricted problem of three bodies, 
the choice of the advances in the orbit-to-orbit progression 
will be guided by the desire of recording all interesting events 
generated in the series within the time period each orbit is 
followed. 


Events of the restricted problem that serve to corro- 
borate the comparativeness of orbits of the two problems 
studied are the occurrences of periodic orbits that can be 
identified as counterparts of Kepler orbits. These three-body 
orbits will be assigned an identification in terms of n*-values 
that are identical to the values characterizing the referenced 
Kepler orbit. The n*-value, of course, will not have the 
physical meaning of relating the mean angular motion as the 
n-value has, though in most cases it may be taken as a fairly 
good approximation. 


However, the series of the restricted problem exhibits 
many interesting events that are peculiar to itself and are 
not to be found in the Kepler series. This is responsible for 
a large number of orbit drawings to be shown for which no 
Kepler orbit counterparts are furnished. 


For deciding on an appropriate starting point in the 
series, knowledge of practical identity of the two series in the 
hyperbolic range precluded showing orbits in this range. In 
fact, this practical identity is preserved for orbits in the 
elliptical range, if they are still close enough to the para- 
bolic boundary. 


Thus, this area provides a natural starting point. The 
first three orbit pairs depicted (Figures 23-28) demonstrate the 
fact of near-identical behavior in the two series. The n-values 
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Yr *- 2.3363 
P/2 = 6.25 


Figure 23 Figure 24 Figure 25 

are -1/4 , -1/3, and -1/2, and correspondingly the n*-values. 

From these examples it is also clear what shape and topology 
the orbits of smaller (absolute) n-values will have (as -1/5, 
-1/6, etc.) approaching ever closer the parabolic orbit. 

The fourth orbit pair illustrated (Figures 29, 32) is that 
of a circular orbit in the Kepler series and its near-circular 
counterpart in the upper series. 

Subsequently, a first example of orbits is shown where 
the Kepler orbit cannot be paired to a periodic orbit of the 

KEPLER ORBITS IN ROTATING COORDINATES 


1 L_ 1 1 1 1 = —2.4215 

n **' ,/4 P/2 = 9.41 

Y R = -2.4657 

P/2 = 12.55 
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MEAN ANGULAR MOTION n»-l/4 
INITIAL VELOCITY Y R »-2.4677 
HALF -PERIOD P/2-12.57 
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n*-l/3 
Y r = -2.4234 
P/2 =9.42 


n =d-l/2 
Yr =-2.3386 
P/2 = 6.28 


ORBITS START AT X R -1.50 


Figure 26 


Figure 27 


Figure 28 
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(CIRCLE - COORDINATED) 
Yr =-2.3138 
P/2 = 2.02 


Figure 29 



(BEST APPROACH TO PERIODIC) 
* R =-2.2455 
REF. TIME : 9.21 

Figure 30 



n * = -I 
'Y'r =-2.1005 
P/2 = 3.21 

Figure 31 


top series (Figures 30 and 33) . The best orbit in the three- 
body . series to match the Kepler orbit of n = -2/3 is that shown 
for Y r = -2.2455, which at the comparable axis-crossing fails by 
0.3 degree of reaching orthogonality. 


The lack of agreement of this orbit with the Kepler orbit 
of n = -2/3 is explainable by the nearness at which this orbit 
passes at its second X R -axis crossing by the Moon — located at 
X R = 0.9875. In contrast to this, the next E-M orbit (i.e., 
that for n* = -1 on Figure 31) has its second X R -axis crossing 
sufficiently away from the masses that its behavior is quite 
similar to its Kepler counterpart (Figure 34) . 


The fact that the two orbits on Figures 30 and 31 pass by 
the Moon on opposite sides when they cross the X R -axis the 
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CIRCULAR ORBIT (NEGATIVE) 
Yr =-2.3165 
P/2 = 2.03 


Figure 32 



n =-2/3 
Y r =-2.2551 
P/2 = 9.42 

Figure 33 
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n =- 1 
Y r =-2.0773 
P/2 = 3.14 

Figure 34 
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REFERENCE ORBIT 
* R =-2.22020 

FLY BY TIME (MOON): 3.720 
FLY BY TIME ( EARTH) : 4.132 



COLLISION ORBIT 
* R =-2.22012 

FLY BY TIME (MOON): 3.720 
COLLISION TIME (EARTH) : 4.119 


Figure 35 



Figure 36 


REFERENCE ORBIT 
* R =-2.2200 

FLY BY TIME (M00N):3.7I9 
FLY BY TIME (EARTH):4.III 
Figure 37 


second time points to the existence of a lunar collision for an 
orbit with intermediate initial velocity. In an attempt to 
isolate this occurrence, one comes across the interesting case 
of two collision orbits being close together, one colliding 
with the Earth, the other colliding with the Moon. These events 
are delineated in a sequence of six orbits in Figures 35-40. 

For the first three of the orbits (Figures 35-37) , attention is 
drawn mainly to the orbit behavior near Earth. The classical 
sequence of events near a collision orbit is developed insofar 
as the second orbit in the group depicts the exact collision, 
while the first and third orbits show close flyby orbits in 
directions around the Earth opposite to each other. For the mean 
angular motion "in the large," however, it is worth noting that 
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REFERENCE ORBIT 
* R =-2.2192 

FLY BY TIME (MOON). 3.715 



COLLISION ORBIT (MOON) 
Y r = - 2.2190 
COLLISION TIME: 3.715 


Figure 38 



Figure 39 


REFERENCE ORBIT 
Y r = - 2.2188 

FLY BY TIME (MOON): 3.714 
Figure 40 


all orbits finally again assume the retrograde motion which they 
originally followed. 

While for the three orbits shown in Figures 35-37, the 
motion around the Moon shows little change, the three orbits 
shown in Figures 38-40 then develop the events connected with 
the lunar collision. The orbit shown in Figure 39 represents 
the exact collision with the Moon. It is flanked by close flyby 
orbits, doubling the lunar mass point in directions opposite 
to each other, as seen by an observer on the Moon. Again, after 
leaving the neighborhood of the mass, all orbits assume the 
retrograde sense of mean angular motion, which they originally 
followed. It is of interest to observe that all six . orbits 
are initiated by velocities that lie in a space of AY R = 0.0014. 


22 



At this point. Figure 41, which is an annotated version of 
the apsidal-velocity diagram (Figure 16, page 12), will serve to 
orient the reader as to the current position in the series. On 
this figure, the range that has been covered up to this moment 
is marked by the lower bracket. 

The bracket above it embraces the range in the series that 
will be entered shortly. As seen on this diagram, the series 
enters the range of its highest n-values. 

The peak is reached with *R T ” 1,5 which is the (only) 
Kepler collision orbit of this series. Past this orbit in the 
development, the mean angular motion, which reaches a value of 
slightly more than |n| = 3/2, reverses its sign and falls off 
in value. 



Figure 41.- Apsidal velocities in the rotating system for orbit families with fixed mean motion 


Figures 42 through 53 depict orbits falling into the n- 
range from -4/3 to +4/3, with the three-body orbits (top of page) 
closely resembling the Kepler orbits (bottom of page) . 
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BEST APPROACH TO PERIODIC 
? R =-1.810 
Tr £F = 9.46 

Figure 42 


COLLISION WITH MOON 
fp = - 1.750 
COLLISION TIME: 5.97 

Figure 43 


Yr = - 1.568 
P/2 = 6.30 

DASHED LOBE CUTS X R -AXIS 
AT 90° 

Figure 44 


In particular, the first and third Kepler orbits (Figures 
45 and 47) show the n-values of -4/3 and -3/2. The orbits on 
Figures 52 and 53 have the same n-values, but of positive sign and 
in reversed order. The orbit on Figure 51 is the collision orbit, 
with its initial velocity value identical, but negative, to its 
position value. The second orbit of this group (i.e.. Figure 46) 
is provided here for reason of comparison with its counterpart of 
the three-body series. 

The orbits of the three-body problem, which are arranged 
(Figures 42, 43, 44, 48, 49, and 50) appropriately above the 
corresponding Kepler orbits, reflect approximately the development 
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Figure 46 


Figure 47 
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n* ■ ±3/2 
=-1.500 

COLLISION TIMES: 2.0; 6.3; 10.6; ETC. 
PERIODIC COLLISION WITH EARTH 

Figure 48 
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REFERENCE - ORBIT 
Yr =-1.35 

T R#f = 6.54 (ON X R -AXIS) 

REF ANGLE: 97° _ 

Figure 49 



*R - - 1.21 
T Rtf * 9.52 

REF ANGLE: 85.7" 
Figure 50 


of the Kepler orbits. The first of the group is a "best match" 
to the Kepler orbit of n = -4/3, while the third is periodic and 
fits well the structure of its Kepler partner. Again, "best 
match" is all that could be found to compare with the last Kepler 
orbits of n = +4/3. The three-body counterpart to the Kepler 
collision case, however, is a fine example of a periodic collision 
orbit, with initial velocity not different from that of the Kepler 
orbit. Its periodicity characteristic is n* = +3/2; thus, this 
orbit combines the characteristics of two of the opposite series. 
Of the remaining two orbits of the upper group, the orbit (with 
Yr = -1.35) following the collision orbit only serves the purpose 
of suggesting the mode of transition between its neighbors. The 
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n = +4/3 
Yr = - 1.151 
P/2 = 9.42 
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Figure 52 


Figure 53 
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REFERENCE ORBIT 
Yr = -1.130 

T Ref= 680 

Figure 54 



PERIODIC ORBIT 
Y r =-1.12884 
P/2 = 6.80 

Figure 55 



COLLISION 
WITH MOON 

Yr = “1.1260 
COLL. TIME: 6.80 


Figure 56 


earlier orbit, however, with Y R = -1.750, exhibits another mode 
of colliding with the Moon. Its neighboring orbits were not ex- 
plored further, though it is obvious that, to its left, before 
the velocity reaches Y R = -1.810, there is again a case of colli- 
sion with the Earth. 


Figures 54-58 represent five orbits of the restricted prob- 
lem. They may be loosely identified as being in the transition 
area between general orbits and lunar satellite orbits. "General 
orbits," in this connotation, can be defined as orbits that have 
a non-zero mean angular motion with respect to the heavier mass 
as reference point, and lunar satellite orbits are orbits of 
zero angular motion with respect to the heavier mass, but non- 
zero angular motion with respect to the smaller mass as reference 
point, with all motion observed in the rotating system. 

This group of five orbits (Figures 54-58) shows that there 
is no sharp dividing line between these two classes. The second 
and fifth orbits of the group are lunar satellite orbits, while 
the fourth is meeting the characterization of a general orbit. 

All three are periodic and therefore known in their full courses. 
The very first of this group is not identified in this respect; 
its showing serves the purpose of giving aid in recognizing the 
transition between its neighbors. The reason for including the 
middle orbit of the group (with Y R = -1.1260) is obvious, viz., 
that this is another collision orbit with the Moon. 

The two Kepler orbits shown in Figures 59 and 60 illustrate 
the trend that is characteristic for the current velocity level; 
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PERIODIC ORBIT 
* R □ -1.0997 
P/2 = 14.95 


PERIODIC ORBIT 
? R □ -1.0604 
P/2 = 12.21 


Figure 57 Figure 58 

i.e., that the angular spacing between successive loops decreases 
with progression in the series. This fact supports materially 
the formation of lunar satellite orbits in the restricted problem 
of three bodies. 


The region of lunar satellite orbits is the subject of 
interest for the next six orbits shown in Figures 61-64, 67, 68. 
First, however, the development that is in progress in the com- 
parable region of the Kepler series may be studied. This develop- 
ment is represented in Figures 65, 66, 69, and 70. The decrease 
in angular spacing of successive loops, observed on the former 
diagram of the Kepler series, has progressed further with the 
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REFERENCE ORBIT 
Yr = -1.100 
T R»f * 6 80 

Figure 59 



REFERENCE ORBIT 
Yr = -1.060 

T Ref = 680 

Figure 60 
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LOBE-HALF TIME (MEAN):3.03 

(PERIODICITY NOT ESTABLISHED) 

Figure 61 Figure 62 

orbit in Figure 65. The n-value here is +16/15, which indicates 
that there are 16 loops in the orbit. An infinite number of 
orbits of n-values constructed as ratios p/ (p - 1) with increasing 
integers p will approach a limiting orbit whose n-value is unity. 

A comparison of this limiting orbit, however, with the orbit of 
n = 1 that is plotted in Figure 66, must register discrepancies 
in geometry as well as orbital period. A more appropriate embed- 
ding of the pictured single-loop orbit into its neighboring orbits 
would perhaps be obtained if in each of the neighboring orbits 
just one loop would be considered. Continuity in shape as well 
as in time-spent would prevail, but, of course, periodicity does 
not exist for these neighboring loops. 

The single-loop Kepler orbit of n = 1 in Figure 66 is 
interesting for various reasons. It belongs to the family of 
orbits all of which include the point X R = 1 and Y R = 0 and they 
do not intersect each other. Their mean angular motion with 
respect to the mass, when taken in the rotating system, is zero. 
Orbits of this family have the character of libration orbits, 
with the point X R = 1 and Y R = 0 understood as the libration point 
of the Kepler system in the normalized rotating coordinate system. 
(Of course, any point at unity radius does as well.) 

The correspondence between Kepler orbits and three-body 
orbits in the current energy region is rather weak. A remarkably 
good case of comparability, however, is found in the pair of 
orbits that consists of the Kepler libration orbit and its match 
in the upper series, shown on Figure 64. This latter orbit re- 
produces the geometry, velocity, and orbital period of its Kepler 
counterpart rather well. The similarity between them raises the 
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P/2 = 2.973 


LOBE-HALF TIME (MEAN):3.00 
(PERIODICITY NOT ESTABLISHED) 

Figure 63 Figure 64 

vexing dilemma of association. This orbit is considered to be a 
lunar satellite orbit since it meets the definition, but it is 
recognized that it owes its main characteristics to the gravity 
field of the Earth. 

The Kepler orbits depicted in Figures 6 9 and 70 reveal that 
the trend in the orbital progression reverses now, manifested in 
the increase of loop spacing from orbit to orbit. Also, loops 
are now falling behind, over time, instead of moving ahead, as 
observed before. Thirdly, while loops before were connected to 
each other at the "inner envelope," now the loops are beginning 
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n =+16/15 n *+ 1 

=-0.96207 =-0.922650 

P/2 = 47.12 P/2 n 7 r 

LOBE HALF-TIME n2.95 
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P/2 = 20.1 P/2 = 8.39 

MEAN LOBE -HALF TIME = 2.87 MEAN LOBE - HALF TIME = 2.80 

Figure 67 Figure 68 

at the outer envelope and, after the loop-formation, are return- 
ing there to join the succeeding loop. 

The orbits of the upper group show similarity with the Kep- 
ler orbits in their initial orbit geometry while otherwise they 
develop rather independently. The progressively increasing com- 
pactness of the first three orbits, however, parallels the trend 
that is predominant in the Kepler group (Figures 61-63) . 

No special effort was made to isolate periodic orbits here. 
To the right of the single loop orbit, two periodic orbits were 
encountered which are reproduced in Figures 67 and 68. The period 
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n =+23/24 
Y r =-0.89891 
P/2 = 75.4 

LOBE - HALF TIME = 3.28 
Figure 69 



n =+15/16 
Y r =-0.08726 
P/2 □ 50.27 

LOBE -HALF TIME = 3.35 
Figure 70 
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PERIODIC ORBIT 
=-0.879860 
P/2 = 25.98 


PERIODIC ORBIT 
'Y'r= -0.875970 
P/2 = 21.71 


Figure 71 


Figure 72 


PERIODIC ORBIT 
Yr= -0.874621 
P/2 = 18.25 
Figure 73 


covers seven loops on the earlier orbit and three loops on the 
latter orbit. 


Next in the three-body series is a triple of closely 
neighbored orbits (Figures 71-73) that illustrate one possible 
mode of "losing a loop" in an orbital development. The loop to 
be lost is the one shown in dashed lines in Figure 71. This 
loop travels along the orbit as one progresses in the series. 

Figure 72 shows this loop placed at the negative Xj^-axis. In 
Figure 73 the loop is combined with the loops present in the 
lunar area, without adding to their number. 

Going back to Figure 71, one finds that during this pro- 
gression actually two loops vanish, the second one being the 
mirror image of the first one. This latter one vanishes exactly 
with attainment of the second orbit of the triple. 

This group of orbits is noteworthy for a second feature as 
well, viz., the partition of each orbit into a time span of lunar 
captivity with the remaining time spent in a more general orbit. 

It is interesting to see that the breakaway is in most cases 
accomplished with a close flyby of the Moon. The reverse pro- 
cedure, the capture by the Moon, is accomplished in a similar 
manner . 

More examples of motion alternating between two modes of 
flight are exhibited in the four orbits shown in Figures 74-77. 

In fact, the motion here alternates between a lunar satellite phase 
and an Earth satellite phase. The latter is understood to have 
zero mean angular motion with respect to the Moon as seen in the 
rotating system. The Earth satellite mode can well be identified 
when comparison is made with the Kepler orbits shown in Figures 11 
and 12. So, for example, the dashed part of the orbit of Figure 
74 shows the characteristics of the Kepler orbit of n = 3 , as 
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Figure 74 Figure 75 

depicted in Figure 11, and the Earth 
satellite portion of the orbit in 
Figure 75 corresponds well to the 
Kepler orbit with n = 4/3, as depic- 
ted in Figure 12 . The next in line 
(Figure 76) resembles that of n = 3/2, 
as seen in Figure 11, and the last of 
the four orbits (Figure 77) suggests 
the characterization by n = 5/3 for 
its Earth satellite portion. 

On all four of these orbits the 
lunar close approach, which triggers 
the transition from mode to mode, is 
quite pronounced. 


Figure 76 



PERIODIC ORBIT 
Y R = - 0.855166 
P/2- 13.33 

Figure 77 


The collision orbit, out of which the first of the four 
orbits is developed, is not shown here. The remaining three 
orbits progress toward another collision orbit. Its image is 
traced on Figure 81. 


Since this lunar collision orbit is the last one of its 
nature, recorded in this report, the cases will be reviewed here. 
If, in the review, those orbits that reach collision after first 
completing a full loop (or more) around the Moon, as e.g., would 
be the case for an orbit between Y R = -.8746 and Y R = -.8700, 
are excluded here, four orbits would remain. The three collision 
orbits shown earlier are replotted here in Figures 78, 79, and 80. 
The velocity direction at which the orbits impact the Moon, 
especially the consistent angular progression of it from case to 
case, is noteworthy. Also worth noting is the relative timing of 
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COLLISION ORBIT (MOON) 
* R =-2.2190 
COLLISION TIME: 3.715 


COLLISION WITH MOON 
* R =-1.750 
COLLISION TIME: 5.97 



COLLISION 
WITH MOON 

'Y'r = -1.1260 
COLL. TIME: 6.80 


Figure 79 



Yr = -0.8489 
COLLISION TIME : 4.18 


Figure 80 Figure 81 

the collision events with the last one shown not being the earli- 
est in time, as might perhaps be suggested by the short path length. 


If the significance of a collision orbit is measured by 
the number and variety of events that are brought about in the 
flyby orbits in the neighborhood of the collision orbit, then 
the last collision orbit is ranking highest in this respect among 
the four listed. As the subsequent figures (82-99) will prove, 
many types of orbits are packed into the space generated by the 
small increase of initial velocity of less than Av = 0.001. 


Figure 82 portrays a periodic orbit that, after flyby of 
the Moon, makes an extended excursion into the field ahead of 
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PERIODIC ORBIT 
=-0.848798 
P/2 = 29.70 


Figure 82 



REFERENCE ORBIT 
f R =-0.84856 
T Ref = 23.4 

Figure 83 



REFERENCE ORBIT 
Y r =-0.84850 
T Ref =15.84 

ANGLE AT T Raf = 88.7° 
Figure 84 


the Moon. From there it returns to the Moon, crossing the X R 
axis orthogonally between the two masses. The insert in Figure 
82 delineates the flyby arcs in enlarged scale. In its full 
period, this orbit makes only one revolution about the Moon, with 
the sense of the revolving motion being direct. The (normalized) 
time period is nearly 60 units which corresponds to about 9 
months in the real Earth-Moon system. 


Note that the loops are numbered according to their sequence 
in time. So are the loops of the next six orbits. The purpose is 
to draw attention to the relative position of these loops and the 
formation through which they are going over the course of the 
subsequent orbits. 


Figure 83 demonstrates the beginning of the formation with 
the loops extending forward in sequence of their numbers. On the 
subsequent orbit (Figure 84) loops are spaced far apart, the 
fourth loop extending into the lunar area. Here, also, it may be 
recognized that with little effort this orbit can be modified to 
let the crossing at T = T Re f be exactly orthogonal whereby a 
periodic orbit is obtained that encloses both Earth and Moon. 


In Figure 85, the number 2 loop is coming close to the 
negative X R -axis. In Figure 86, this loop is straddling the 
axis, cutting it at 89.4 degrees. A minute correction in the 
initial velocity will result in an orthogonal crossing at this 
place. (It is worth observing that the correction is concerned 
with the fifth significant digit!) 
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REFERENCE ORBIT 
Y r =-0.84847 
T R „= 19.25 

ANGLE AT T„„ = 87.0'’ 
Figure 85 



REFERENCE ORBIT 
* R =-0.84845 
T| = 12.9 WITH < =94.6° 
T 2 b 23.0 WITH < =89.4° 

Figure 86 



REFERENCE ORBIT 
Yr =-0.84840 
T R ,f= 19.4 

ANGLE AT Tr„ = 103.8° 
Figure 87 


In progressing to the sequent orbit (Figure 87) a number 
of new loops is coming into existence (six may be counted) . 

Here again the possibility of isolating a periodic orbit is at 
hand, since now loop number 5 is straddling the negative X^-axis . 


Other chances for formation of periodic orbits are found 
between two straddling positions by the fact that the arcs link- 
ing two loops may furnish the right angle. Such a situation is 
approached, for example, in Figure 88 where the currently measured 
angle of 89.6 degrees indicates the proximity of a periodic or- 
bit . 


To introduce the orbit that now follows in the progression, 
another look at the 'orbit on Figure 82 is suggested. That orbit 
is the only one so far in the current group that can be con- 
sidered to be of a lunar satellite mode. The orbit in Figure 89 
now falls into the same classification. The topological charac- 
ter of this orbit is identical to that of Figure 82 insofar as 
it actually forms one closed loop around the Moon. This can be 
verified by means of the blown-up portion of the insert. The 
period of the present orbit, however, is much shorter than that 
of the former one. The character of the orbit changes drastically 
if the initial velocity is raised by 1/1000 of a percent as is 
borne out when passing from the last orbit to that of Figure 90. 
The new one encompasses the Earth in a wide near-circular arc 
after which it returns to the lunar area. This orbit is not 
developed to its exact periodic shape, but the angle of 88.6 
degrees at the axis crossing at T^ e f points toward the nearness 
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t Ref = 15.0 WITH <£=88.6“ 


Figure 88 


Figure 89 


Figure 90 


of this event. It is interesting to determine its topological 
classification (in the completed periodic shape) . After elimina- 
ting all loops that do not encompass a masspoint, it will be 
found that the topological structure is that of a lemniscate. 

Between the two orbits of Figures 90 and 91, there is an- 
other orbit that collides with the Moon. The collision occurs 
when the Moon is approached the second time. The collision is 
evidenced by the fact that in Figure 90 the orbit arc that makes 
the second pass by the Moon crosses the X^-axis "at the left" of 
the Moon, i.e., between the Moon and the Earth, while in Figure 
91 this crossing takes place "at the right" of the Moon. 

The little loop which has formed around the Moon on the 
last orbit (Figure 91) will grow and decisively determine the 
development of the following orbits. The process shown in the 
next two of these orbits (Figures 92 and 93) is self-explanatory. 
In Figure 93, the loop referred to draws parallel to the first 
loop, and, as shown in Figure 94, these two loops have combined. 
This then is an orbit that does not enclose either of the masses. 
It is associated with the translunar libration point, designated 
Li in reference 1. Past this L^-libration orbit, the second 
loop shows up again, now being outside and below the first loop. 
This is shown by the dashed curve of Figure 95. This loop 
"rolls away," while the progression in the series goes on. 
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REFERENCE ORBIT 
Y r =-0.84830 


REFERENCE ORBIT 


REFERENCE ORBIT 
Y r =-0.8479 


Figure 91 


Figure 92 

EARTH-MOON FIELD ORBITS 


Figure 93 


LIBRATION POINT 



PERIODIC (LIBRATION ) ORBIT 
Yr =-0.847844 
P/2 = 4.1952 


REFERENCE ORBIT 
Yr =-0.847800 


PERIODIC ORBIT 
Y r " -0.847327 
P/2 = 20.48 


Figure 95 


F igure 96 


The general development is rather intelligible at this 
point in the series. A number of periodic orbits is met 
(Figures 96-99) which, on their Earth satellite portion, follow 
a pattern that is also encountered for Kepler orbits in the 
current velocity range. For comparison, the behavior of Kepler 
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PERIODIC ORBIT 
Y r =-0.847252 
P/2 = 17.31 



PERIODIC ORBIT 
Yr o -0.847091 
P/2 = 14.11 



PERIODIC ORBIT 
Y r =-0.846990 
P/2 = 20.36 


Figure 97 Figure 98 Figure 99 

orbits referred to is illustrated in Figures 100 and 101 below 
and 104 and 105 at the bottom of the next page. These figures 
demonstrate that loops are getting smaller and drifting apart. 


The four orbits plotted in Figures 102, 103, 106, and 107 
represent a group on its own, of which the principle of develop- 
ment can best be recognized by studying them in reverse direction, 
starting at Figure 107. This orbit shows a small loop in the 
lunar area, with the orbit encompassing the Earth once within 
its period. With a move away from this orbit to those on its 
left (i.e. , earlier orbits), the small loop is increasing and 
also approaching the Moon. The "energy-gain" in this successively 
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REFERENCE ORBIT 
Yr *-0.8430 

T Rtf - 9-334 


PERIODIC ORBIT 
Y r » -0.830738 
P/2 « 16.12 


Figure 102 Figure 103 

closer flyby mode causes the orbits to increase their overall 
diameters, by which then also the number of revolutions per 
period about the Earth increases. The orbit of Figure 102, which 
is the last shown of this development, appears to be an escape 
orbit. (For later reference, this group will be called the 
"pigtail" group.) 

Comparing the orbit of Figure 105 here with the orbit of 
Figure 107 indicates that the small loop in the three-body-orbit 
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n =+5/6 
Y r —0.830833 
P/2 « 18.84 

Figure 104 


n =+4/5 
Yr =-0.813298 
P/2 = 15.71 

Figure 105 
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PERIODIC ORBIT 
Y r =-0.827640 
P/2 * 12.88 


PERIODIC ORBIT 
Y r =-0.815682 
P/2 * 9.460 


Figure 106 


Figure 107 


is a "Kepler-characteristic," but that the effect exerted by 
the Moon on the loop during flyby is strong enough to change the 
flight pattern following the loop. 

While in the velocity bracket between Y R = -1.2 and 
Y r = -0.8, which has just been covered, the two series compared 
showed a rather high disparity — except at the very beginning 
of each orbit and excepting the singular case of good similarity 
found between the Kepler libration orbit and the single-loop 
lunar satellite orbit — the orbits following from hereon exhibit 
a gradual restoration of comparability. Similarity will first 
be found in the topological structure, later also in the geometry. 

Before discussing Figures 108 and 109, the progression of 
events may be picked up first at the Kepler series as shown in 
Figures 110-112. The three orbits of Figures 110-112 effectively 
portray the developments on both sides of a cusping orbit, with 
loops seen before it and smooth indentation past it. Notice 
that the n-value of a cusping orbit, in general, is incidental, 
like that of a circular orbit. It is a function of the initial 
distance from the origin. 

Now, with respect to the three-body orbits shown in Figures 
108 and 109, coordinated to the cusping Kepler orbit, one notices 
that the development follows the same line, but that one cannot 
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CUSPING ORBIT 
Y ft =□ 0.7600 

TIMES: T, = 4.02; T 2 = 12.25; T 3 =^21.50 
Figure 108 



CUSPING ORBIT 
Y r = - 0.7500 

T, = 3.76; T 2 = 12.77; T 3 = D 21. 94 
Figure 109 


speak of a single cusping orbit. Rather, each of the loops that 
occurs over time on the three-body orbits observes its own 
"maturing process", thus spreading the occurrence of cusps over 
a range in the progression. 
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n »+3/4 
Yr *-0.78749 
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Figure 110 



CUSPING ORBIT 

Y r *-0.75801 

T Cu *p * 4.54; 13.62; 22.70 


Figure 111 



n * + 2/3 

Yr - -0.74489 
P/2 - 9.425 


Figure 112 
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n* = + 5/8 
fp = -.73839 
P/2 = 27.0 

Figure 113 


H* o + 3/5 
Yp = -.72640 
P/2 = 17.4 

Figure 1 1 4 


In Figures 113 through 120, orbit correlations and geo- 
metrical similarities are noted. In Figures 113-116, the 
comparability of the two series reaches the point where orbits 
can be correlated with respect to their n (n*) -values . Geo- 
metrical similarity is accomplished favorably in Figures 117- 
120. The first orbits in these figures are circular in the 
Kepler series and near-circular in the upper series. The second 
orbit (Figure 120) is of value n = +1/2 in the Kepler series; 
the counterpart in the upper series fails by 0.8 degree of 
paralleling the periodicity. 
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n =+5/8 
? R = -.72390 
P/2 □ 25.1 

Figure 115 



n =+3/5 
? R =-.71136 
P/2 = 15.7 

Figure 116 
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FROM: POSITIVE CIRCULAR 


BEST APPROACH TO PERIODIC 


* R =-.71854 
P/2 = 8.434 


* R a -.6854 
Tref- 7.20 


Figure 117 
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From this point, there is a good one-to-one correspondence 
of Kepler orbits to those of the restricted problem of three 
bodies, at least as long as orbits are compared whose values 
of n (or n*) have low integers in the numerator and denominator. 
Test cases shown concern n = +5/11, +4/9, +3/7, and +2/5, all of 
which are plotted in Figures 124, 125, 126, and 130, with cor- 
responding n*-cases in Figures 121, 122, 123, and 127. 

The arithmetic progression in n-values would ask for the 
two cases of n = 1/3 and n = 0, the last one being the parabolic 
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CIRCULAR (POSITIVE) 
= -.68350 
P/2 ■ 6.89 


n = + 1/2 
= - .66133 
P/2 = 6.28 


Figure 119 
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PERIODIC ORBIT 
n* = +5/11 

= -0.640597 
P/2 = 34.77 


PERIODIC ORBIT 
n* = +4/9 
f R = -0.632473 
P/2 = 28.48 

Figure 122 


PERIODIC ORBIT 
n* b +3/7 
Y r = -0.624343 
P/2 = 22.26 

Figure 123 


escape orbit. The parabolic orbit was omitted because the illus- 
tration of orbit loops continuously spiralling-out is well known. 

Instead, the orbit of n = 1/3 is depicted together with that 
of n = 1/4 (Figures 131 and 132) . This serves to show the begin- 
ning development of the infinite sub-series of orbits, all associ- 
ated with n-values of the type 1/1 (with increasing integers £) 
that converges toward the parabolic orbit. A parallel develop- 
ment holds for the E-M orbits , of which the cases of n* = 1/3 
and n* = 1/4 are shown in Figures 128 and 129. 

KEPLER ORBITS 


Yr = 
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Yr * 
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n = + 1/3 
* R = -0.57665 
P/2 * 9.42 



n * + 2/5 
? R - -0.610927 
P/2* 15.71 

Figure 130 


n = -t- 1/4 
Yr = -0.53228 
P/2 = 12.57 


Figure 131 


Figure 132 
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CUSPING ORBIT 
= 0.00 


REFERENCE ORBIT 
* R =+0.300 


Figure 133 Figure 134 

The two final orbits that are presented for either series 
(Figures 133-136) serve to illustrate events and developments 
encountered in the hyperbolic range. The orbit of Figure 135 
represents the only other cusping orbit there is in the Kepler 
series. The cusp is at the beginning of the orbit where the 
velocity coordinates are X R = Yr = 0. For better visibility of 
the cusp, some negative time is shown on the orbit (dashed curve) . 
The series of the restricted problem of three bodies contains the 
exact counterpart to this orbit. 


The orbits of Figures 134 and 136 exemplify the development 
the two series take after passing the previously mentioned cusp- 
ing orbits. The cusps are succeeded by loops that are growing in 
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CUSPING ORBIT 


fR =0.00 


tp = + 0.300 


Figure 135 


Figure 136 
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size with progression of the series. The symmetry points of the 
loops are at the points of orbit initiation. The initial velocity 
is positive from now on. However, consistently, the angular motion 
on all these orbits will eventually reverse to the retrograde 
direction and remain so. 


National Aeronautics and Space Administration 
Electronics Research Center 

Cambridge, Massachusetts, August 1968 
129-04-04-06 
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APPENDIX 


DIFFERENTIAL EQUATIONS AND 
METHODS OF REGULARIZATION AND INTEGRATION 
EMPLOYED IN THE COMPUTATIONAL PROGRAM 


1 .- REGULARIZATION OF THE RESTRICTED PROBLEM OF THREE BODIES 

Part 1 of the Appendix considers the equations of motion of 
the planar restricted problem of three bodies . The Newtonian 
equations of motion are regularized to eliminate singularities 
and to allow the analytic continuation of solutions passing 
through the singular points. The regularizing transformation and 
its development is due to Arenstorf (ref. 1 ). 

Consider the motion of three bodies in a two-dimensional 
rotating system. Assume the only forces in inertial space are 
inverse square attractive forces acting between the bodies. 

Assume masses m]_ and mj are in circular orbits about their common 
center of mass and mass m3 is sufficiently small that it does not 
affect the motion of m]_ or m2. 

The planar equations of motion are given in a coordinate 
system with origin at the center of mass of m^ and m2 / and rota- 
ting with respect to inertial space with angular velocity such 
that m^ and m2 remain fixed on the rotating X R axis (Figure A-l) . 
C-| and C2 (C^> 0 , C2 < 0 ) are constant vectors on the Xp-axis and 
locate mi and m2 with respect to their center of mass. The limits 
on Ci and C2 are discussed in Part 2 of this Appendix. r locates 
m3 with respect to the origin; ri locates m3 with respect to mi; 
and r_2 locates m3 with respect to m2 . 
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The equations of motion of m 3 in the above system are 


f + 2w x r 


2 

oj r 



( 1 ) 


with w, the angular velocity of the coordinate system, directed 
normal to the plane of motion, and the magnitude given by 


to 


/G (m^ + m 2 ) 


< c i - c 2 > 


(2) 


2 

The regularization transforms the 1/r singularities so 
that the new system is integrable (in a computational sense) and 
continuous at these points . 

The normalization of the equations of motion could be 
carried out in the usual way. The mass unit is chosen such that 
the sum of masses m^ and m 2 is unity. The unit of length is 
chosen such that (C^ - C 2 ) is unity and the unit of time is 
chosen such that G is unity. With these definitions, the cir- 
cular angular velocity to is also unity. 

However, the regularization carried through below is done 
with the non-normalized equations of motion. 

Transform the Cartesian coordinate system Xr,Yr to a com- 
plex coordinate system in which the complex position vector z is 
given by z = Xr + iYR. For simplicity of notation, denote a 
complex function by the subscript 1 and its complex conjugate by 
the subscript 2. Introduce a complex u plane such that u maps 
onto z by 


z k = f(u k ) (k = 1,2) (3) 

where f (u) is real when u is real. This implies that f (u k ) = 
fk(u) (k = 1,2). In terms of u and f (u) , a Lagrangian for the 
restricted problem of three bodies is: 
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L ^ u k'^k^ = ^ f i^l + iwf i) ( f 2^2 " ioof 2^ 


2Gm n 


2Gm. 


V (f 1 C l^ ^ f 2 C l^ V (f l C 2 ^ ^ f 2 -C 2^ 


( 4 ) 


Since it can be shown that the Lagrangian satisfies Euler's 
equation in terms of the u^'s and u^'s, a Hamiltonian can be ob- 
tained for the system from 


H * = S “k v k ‘ L(u k' i k ) 

k 


where v^., the conjugate momenta, are defined by 


v 


k 


9L 

9u k 


(6) 


H* remains constant and can be identified as twice the Jacobi 
constant . 

Let us transform the time t to a parameter s such that 


s 



1 

g (u 1 )g (u 2 ) 


dt 


( 7 ) 


where g(u) is analytic and real when u is real. g (u) is to be 
chosen such that (a) s is monotonically increasing with t, 

(b) the integral exists everywhere, and (c) s approaches a finite 
limit as t approaches a time of collision. 

A new Hamiltonian H is now defined as 


H = g 1 g 2 (H* - h*) (8) 
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for which the canonical differential equations of u^ and in 
the s domain are given by 


, , x 3H 

u k (s) ” av^TiT 
v k (s > “ - au^TsT 


( 9 ) 


and h* is the value of H* defined at t = s = 0 . 
H becomes 


H = g x g 2 


" v l v 2 . . / v 2 f 2 

+ ICjjI 

U-i 
1— 1 
> 

1 

f ' f ' \ f • 

_l r 2 \ r 2 

- H 

m 


m., 


m. 


- 2G 


+ 


- h* 


V( f l“ C l) ( f 2 -C 1 ) V (f l" C 2 ) (VV 


( 10 ) 


and it can be seen from Eq. (8) that H = 0 for all s. 

f (u) and g (u) are chosen such that H contains no trans- 
cendental or algebraic functions and no singularities which 
cause difficulties in the integration of Eq. (9) . 

Make 


f (u) - C x = s 2 (u) 
f (u) - C 2 = r 2 (u) 


( 11 ) 

(12) 


where r(u) and s (u) are rational functions with zeros of first 
order and common poles. A suitable choice for r(u) and s (u) is 

r(u) - iV c i ' c 2 ( u + ir) 

(13) 

s (u) V^l - c 2 ( u - k) 
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In the u-plane, mass is located at u = ±1 and mass m^ is 
located at u = ±i. Eq. (3) is then a 4-to-l mapping or the u- 
plane into the z-plane since ±u and ±u“-^- are solutions for the 
same z . 


Prom Eq. (12), whenever s (u) or r(u) are zero, f'(u) must 
also be zero. Therefore it is required that g(u) has first- 
order zeros corresponding to the zeros of f’ (u) , r(u), and s (u) . 
Furthermore, g(u) must satisfy conditions (a), (b) , and (c) of 

Eq. (7). 

A choice for g(u) is 


g(u) = uf'(u) = r(u)s(u) (14) 


Finally, H in its simplest form is 

H = U 1 V 1 U 2 V 2 + iu)(u 2 V 2 f 2 g l " u l v l f l g 2 ) 
" 8G(m l r l r 2 + m^s,,) - g^h* 


(15) 


Once the initial conditions in the z-plane have been trans- 
formed to the u-plane, Eq. (9) may be integrated together with 


t 



g (u x )g (u 2 ) 


ds 


(16) 


and the result transformed back to the z-plane. 

The only singularities occur when u = 0,°°; but from Eqs . 

(12) and (13) , these correspond to z = 00 and hence t = 00 . The 
singularities therefore present no problem in the actual computa- 
tion of a trajectory. 

For reasons to be discussed in the following section, the 
equations were integrated by power series. Integration of Eqs. 

(9) and (16) by series in u and v is impractical, thus motivating 
the development and integration of a system of equations explic- 
itly independent of u and v. 
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Define the parameters 


a 1 = r (u 1 ) 

&1 = s (u x ) 

Y 1 = 2f(u 1 ) - (C 1 + C 2 ) 

6^ = a-^8-^ (17) 

e l = u l v l 

= e 2 ~ iw6 2 (Yj + + C 2 ) 

n l = a 2 6 l 


Differentiating with respect to s (' denotes d/ds) , making 
use of Eq. (9) , and appending Eq. (16) results in 


a = 


3 i 


Y- 


6 i 


p 1 

e l 


4 




3 1^1 


a l U l 


46 l y l 


= YtP 


1 M 1 


= 10) 


La)[4e 1 (S 1 6 2 - e 2 Y l^ Y 2 + c i + C 2^ 

+ 8G(m^n-j_ + m 2 r * 2 ^ + 4 Y^<5 2 h* 

e' 2 - iw[Y 2 (Y 1 + C 1 + C 2 )y 2 + 46 1 6 2 u 1 ] 

3 1 3 2^2 + V^l 


(18) 


t' = 


46^2 . 
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It should be noted that each equation, except for t, 
represents two equations : one equation for the real part and 

one equation for the imaginary part. Furthermore, closed form 
recursion formulas for the coefficients of the series may be 
developed, thus eliminating numerical differentiation. 

The initial conditions for Eqs . (18) are developed from 

Eq. (17) . Eq. (13) is solved for one of the four values of u in 
terms of the initial conditions defined in the z plane, v is 
determined from Eq. (6) for regions "not too close" to a mass 
from 


v 


1 



ito(y 0 + C n 


C 2 ) 


] 


(19) 


For close approaches to a mass Eq. (19) is not useful 
since z -*■ °°. Under this condition Eq. (15) suggests that the 
initial conditions be derived from u, h* , and one component of 
v. Since H = 0, the other component of v can be easily found. 


2.- REDUCTION OF THE REGULARIZED THREE-BODY PROBLEM 
TO THE REGULARIZED TWO-BODY PROBLEM. 

It is possible to reduce the regularized equations for the 
restricted problem of three bodies to the regularized equations 
for the Kepler problem provided that care is taken in the limit- 
ing process. 

Let us assume that m^ = 0. Since and C 2 are measured 
with respect to the center of mass, which now corresponds to m 2 , 
a first assumption would be that both Ci and C 2 were zero. With 
this assumption, however, the regularizing equations would become 
invalid . 

If mi were zero and m 2 were finite, the equation for the 
center of mass would require that C 2 be zero, but it would allow 
Ci to be arbitrary. To preserve the regularized equations, it is 
necessary that Ci be any positive non-zero number. This defini- 
tion accounts for the limits of Ci and C 2 defined in Part 1 of 
the Appendix. 

The Kepler problem, in a three-body sense, can thus be 
thought of as the motion of a third mass about two other masses, 
one of which has a mass value zero and is located at a non-zero 
distance from the center of mass. 
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3.- NUMERICAL INTEGRATION METHOD AND OUTPUT CONTROL 

In problems of computing families of trajectories, some 
thought must be given to the amount and nature of the computed 
output. It is essential that a minimum of computed points dis- 
play all the required information in a form which can be easily 
assimilated. These considerations — and the ones below — 
motivated the tailoring of the equations so that they may be 
integrated by series . 

In a numerical integration, if no control is placed on the 
number of lines of output, each integration step is bounded only 
by the constraints of computational accuracy. By specifying the 
number of lines of output for any trajectory, additional problems 
are encountered. For example, if too many points are computed, 
some would have to be selectively suppressed, while if not enough 
are computed, more points would have to be generated. The prob- 
lem then is to develop a scheme which will print out an entire 
trajectory in a fixed number of lines and relate all the desired 
information. 

One of the first complications arises since the integration 
is done in the s domain. In the t domain, the integration starts 
at t = 0 and terminates at t = T max , where T max is a specified 
terminating time. However, from Eq. (7) one can see that the 
value of s corresponding to T max is not unique but rather a func- 
tion of the trajectory on which it is computed. Hence, T max must 
be retained as the terminating condition and any scheme based on 
s cannot be used without knowledge of the trajectory. Further- 
more, attempting to relate the number of lines to intervals in 
the time domain would require undesirable iterations in the s 
domain in order to find the value of s relating to a particular 
time t. 

It was resolved that the trajectory had to be computed in 
its entirety to obtain sufficient information such that the out- 
put could be generated within a fixed number of lines. The series 
approach was considered to be the most desirable method of inte- 
gration for several reasons. First, the trajectory could be inte- 
grated faster and could use larger steps than with the standard 
Runge-Kutta techniques. Truncation error in the series could be 
handled in an open-loop fashion by adjusting the interval length 
of the truncated series such that the remainder is bounded. 

Second, in an interval, the series gives a functional representa- 
tion of the trajectory, in contrast to the discrete representation 
given by a Runge-Kutta method. Furthermore, in a system of 
coupled differential equations, once the coefficients of the series 
solution of .a single variable are known, that variable may be 
treated independently of the solutions of the other variables of 
the system. 
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A FORTRAN IV program was written in double-precision which 
integrates Eq. (18) by series. It was found that the absolute 
values of the coefficients were divergent, so the series were 
truncated after 15 terms . The applicable interval length for 
the truncated series was determined such that the sum of the 
absolute values of the last three terms of each series was less 
than 10- 10 . 

Since this will not guarantee convergence, the problem of 
accuracy was attacked by a number of numerical comparison tests 
using a Runge-Kutta method. However, in a practical way, the 
accuracy was checked most reliably with the Kepler cases, since 
these, for uniformity of procedures, are calculated by the very 
same computational program that computes the restricted three- 
body cases . 

For the purposes of the output control, the coefficients 
and the respective interval lengths are saved until the entire 
trajectory is computed. 

With respect to the control of the distribution of points 
to be printed out, it is desired to have a scheme that assigns a 
judiciously denser output to events which occur on a smaller time 
scale than to events which occur on a larger time scale, and yet 
be flexible enough to enhance the dominant characteristics of the 
entire trajectory. 

Having computed the trajectory in its entirety, assigning 
an approximately equal number of printout lines to each interval 
results in a spread of printouts that is not only economical but 
also satisfies the requirements of illustration. For it is the 
very process of finding the interval length of the series and 
determining the corresponding interval of time as a function of 
the location on a trajectory (the last equation of Eq. (18)) which 
relates the weighting of the printouts to the dynamical properties 
of events on the trajectory. 


4.- MAPPING GEOMETRY 

In this part of the Appendix the geometry of the mapping of 
the z-plane into the u plane will be described, with special 
attention given to the images of the X R and Y R -axes in the u plane. 
Reference to this will be made in the following section with re- 
gard to the isolation of periodic orbits and collision orbits. 

Insight of the mapping of z into u is gained by investiga- 
ting the equations mapping u into z. Writing u in polar form, 

u = Re 10 (20) 
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the mapping of u into z is given by Eq. (3) 


z = f (u) 

- {Km- 2 ♦ k° s 26 + k + ^ 

+ i( Cl 4 C2 ) (h 2 - 4r)sin 20 


( 21 ) 


Note first that R and R ^ map into the same z. Thus each 
z will give two values of R except when R = R~^ = 1. 

Mapping the line Y R = 0 into u maps the X R -axis into u. 

Y R = 0 for 0 = n(Tr/2) (n = 0,1, 2, 3, 4) or R = 1. A mapping of 
the X R -axis is given in Figure A-2. 



Figure. A-2 
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a) 


For 0 


0 


C 1< X R <C° 

b) For 0 = tt/2 

C 2**K*- 

c) For R = 1 

C 1>X R > C 2 


1SRS» or 1 > R > 0 


1 < R < CO or 1JR>0 

o*0* J 


The symbol — ► means the mapping of z into u. Also the 
limits written to the left of X R map into the limits written to 
the left of R or 0 and limits written to the right of X R map 
into limits written to the right of R and 0 . 

The Y R axis (not shown in Figure A-2) can be mapped from 
X R = 0 by 


cos 



( 22 ) 


For the interval 0 < 0 < tt/2 , several observations can be 

made : 


(1) 

If 

l c ll 

- 1 c 2 1 ' 

Yr maps onto a 45° line. 

(2) 

At 

to 

R = 
the 

0 and R : 
45° line 

= 00 , the mapping of Y R is asymptotic 

• 

(3) 

If 

|Cil 

< 1 c 2 1 ' 

0 < 45° 

(4) 

If 

1 Cl 1 

> 1 c 2 1 ' 

0 

LO 

A 

CD 


Thus the entire X R axis and the entire Y R axis are mapped 
into the first quadrant of u. After mapping the two half-planes 
Y r >. 0, Y R iO into u for 0 1 0 1 tt/2 
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(1) 0<Y r <~ — ► 1 <R< » 

(2) -oo < Y r < 0 — ► 0 <R <1 


it can be seen that the entire z plane is mapped into the first 
quadrant of u. z may be mapped into the other three quadrants 
to complete the one-to-four mapping by rotations and reflections 
of the first quadrant as z maps into -u by a 180° rotation of the 
first quadrant into the third, and z maps into ±u“l by reflecting 
the first and third quadrant about the real u axis. Care must be 
taken in reflecting since the mappings for (1) and (2) are inter- 
changed as shown in the fourth quadrant of Figure A-2. Finally, 
in the u plane, m^ is located at u = ±1 and m 2 is located at 
u = ±i. 

It may be noted that as long as and C 2 do not simul- 
taneously go to zero, the mapping has significance and its general 
properties are unchanged. Therefore the mapping is applicable to 
the Kepler problem as well as the three-body problem, provided 
the limiting process is as described in Part 2 of the Appendix. 

Trajectories in the z plane can now be related to trajec- 
tories in the u plane. 


5.- TECHNIQUE OF ISOLATING PERIODIC ORBITS AND COLLISION ORBITS 

In all cases considered in this report, the initial condi- 
tions are such that the trajectory starts on the X R -axis with 
velocity orthogonal to the X R -axis. If the trajectory has the 
above initial conditions, a necessary and sufficient condition 
that Eq. (1) be periodic is that the trajectory return orthogonal 
to the X R -axis. The problem of isolating periodic trajectories 
can thus be characterized to finding the magnitude of the initial 
velocity such that the trajectory returns orthogonal to the X R - 
axis on the n^ crossing of that axis. The crossing criterion is 
established to distinguish different types of periodic trajec- 
tories . 

The orthogonality condition is satisfied if X R vanishes 
at the n t ^ crossing. From Eq. (19) : 


• 1 [\ 

Z 1 2[ 6 l f 


- ia)6 2 <5 


l( Y l 


C 1 + 



( 23 ) 


The condition that X R vanishes is that the real part of the 
numerator of Eq. (23) vanishes. This condition simplifies to 
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the vanishing of the real part of 6i£2> since y-j_ (the imaginary 
part of y) vanishes on the X R axis. Using this condition not 
only isolates orthogonal crossings, but it also isolates colli- 
sions since 6 (=g (u) from Eq. (17)) vanishes identically at a 

collision. 

Therefore, the vanishing of the real part of 6]_e2 on the 
n th crossing of the X R -axis isolates orthogonal crossings and 
collisions . 

Let <J> equal the real part of 6i£2 the nth crossing of 
the X R -axis. The following discussion is directed to the 
orthogonal crossing, but may be applied to the collision with- 
out loss of generality. 

The iterative process on the initial velocity first 
establishes the sign of the derivative of the initial velocity 
with respect to <J> . The initial velocity is varied such that the 
absolute value of cp is decreasing until a change of signs in 
occurs. If the sign of the derivative defined above changes 
signs prior to cp changing signs, the initial velocity is iterated 
such that the absolute value of 4> is minimized. This process 
isolates the trajectory which best approaches an orthogonal 
crossing and shows that an orthogonal crossing may not exist for 
the neighborhood of the iterated initial velocities. 

If <)> changes signs, it is assumed that an orthogonal 
crossing exists and that an upper and lower bound on the initial 
velocity has been established. The iteration reduces the inter- 
val in such a way that | <p | approaches zero and is terminated 
when the variation of the initial velocity is in the eleventh 
significant digit. 

This then describes the general approach to isolating 
periodic and collision trajectories. The remainder of this 
discussion will be devoted to the methods used to determine a 
crossing of the X R -axis and the method of evaluating <p . 


6.- CROSSING OF THE X R -AXIS 

The criteria for determining a crossing of the X R -axis will 
now be established. In the z plane, a crossing is indicated by 
a change of signs of Y R . Therefore, let changes of signs of a 
parameter or parameters, as yet unspecified, indicate a crossing. 
Second, the parameters to be tested are constrained to those 
integrated in Eq. (18) or a combination thereof. Third, the 
testing should be done only on those points computed in the ini- 
tial integration of the trajectory. These points are then the 
initial points and final points of an interval, the final point 
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being identical to the initial point of the following interval. 
Fourth, the parameters to be tested should also indicate the 
possibility of a collision. And fifth, the tests should be as 
few as possible and the parameters as simple as possible. 

From Eqs. (17) and (21), y would appear a suitable candi- 
date as a testing parameter. y^, as stated previously, is 
directly proportional to X R and changes of signs in y^ would 
indicate the crossings. The value of y r (the real value of y) 
relative to C-^ or C 2 at the crossing would indicate the proximity 
to a collision. 

However, y does not prove to be a very good parameter for 
the following reasons. Consider in the u plane a single interval 
with the initial point in the first quadrant outside the unit 
circle and the final point in the fourth quadrant inside the unit 
circle (points A and B in Figure A-2) . It can be seen that both 
points lie in a mapping of Yr > 0 and hence y^ would not indicate 
a crossing, while, in fact, either two crossings or a collision 
occurred. Furthermore, y^ could not indicate a collision, since, 
in a strict sense, a collision trajectory touches the X R -axis at 
only a single point and does not cross it, except for collisions 
symmetric about the X R axis. Therefore, y is excluded as a 
useful parameter. 

Although the u plane does exhibit the proper number of 
crossings, we have attempted to remain explicity independent of 
u. Utilizing u at this point would have required twice the 
number of equations and thus twice the computation time to arrive 
at the same result as derived by the method presented below. 

If one had attempted to rationalize the choice of parameters 
in Eq. (17) , one would conclude at first that r) was chosen because 
it simplified the differential equations and the recursion formu- 
las. However, in the light of the past discussion, ri has a much 
greater significance which is brought out in the mapping of u, in 
polar form, into the ri plane. 


n = n r + in i 

Ci ” ^2 / 2\ C x - C 2 

= -=-g —\R - + i -±~2 sin 26 (24) 


Note first that n r changes signs whenever the trajectory crosses 
the unit circle (R = 1) in the u plane and, second, that 
changes signs whenever the trajectory passes through a phase 
angle 6 = n(Tr/2) (n = 0,1, 2, 3, 4) in the u plane. From the 
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mapping of the X R axis in the u plane, it can be seen that the 
change of signs of n r or ni corresponds to a crossing of the X R 
axis. Third, and equally significant, the simultaneous vanishing 
of r) r and m corresponds to a collision. 

The mapping in effect takes the X R -axis mapped in the u 

plane onto the ri r axis and the rii axis and maps the four mass 

points in u into the origin in ri . 

The problem of the double crossing encountered by y. is 

now taken care of by the simultaneous changing of signs by n r 

and rii in the same interval. However, additional refinements 
must be made in this type of situation since the change of signs 
only occurs at the endpoints of the interval. For the purposes 
of counting crossings, isolation of the n^-h crossing, or isolation 
of a collision, it is necessary to subdivide the interval until 
either the crossings lie in separate subintervals or until a 
numerical collision is established. A numerical collision is 
defined when the length of n lies within'a circle of radius 10“^ 
and is then treated as a true collision. The crossing of the 
origin of n thus motivates the count of one crossing for a 
collision. 

Once the interval containing the n^* 1 crossing has been 
determined, it is simple to determine by iteration the sub- 
interval length for which the n r or rii series, or both, go to 
zero, thus locating the X R -axis. Knowing this subinterval 
length, <f> can be computed directly and the isolation of periodic 
or collision trajectories can proceed as described in Part 5 of 
this Appendix. 


REFERENCE 

1. Arenstorf, R.F.: New Regularization of the Restricted Problem 
of Three Bodies. Astronomical Journal, Vol. 68, no. 8, 1963. 


NASA- Langley, 1968 


30 C-45 


63 



